Epigenome-wide association study of dietary fatty acid intake.

BACKGROUND: Dietary intake of n-3 polyunsaturated fatty acids (PUFA) may have a protective effect on the development of cardiovascular diseases, diabetes, depression and cancer, while a high intake of n-6 PUFA was often reported to be associated with inflammation-related traits. The effect of PUFAs on health outcomes might be mediated by DNA methylation (DNAm). The aim of our study is to identify the impact of PUFA intake on DNAm in the Cooperative Health Research in the Region of Augsburg (KORA) FF4 cohort and the Leiden Longevity Study (LLS). RESULTS: DNA methylation levels were measured in whole blood from the population-based KORA FF4 study (N = 1354) and LLS (N = 448), using the Illumina MethylationEPIC BeadChip and Illumina HumanMethylation450 array, respectively. We assessed associations between DNAm and intake of eight and four PUFAs in KORA and LLS, respectively. Where possible, results were meta-analyzed. Below the Bonferroni correction threshold (p < 7.17 × 10-8), we identified two differentially methylated positions (DMPs) associated with PUFA intake in the KORA study. The DMP cg19937480, annotated to gene PRDX1, was positively associated with docosahexaenoic acid (DHA) in model 1 (beta: 2.00 × 10-5, 95%CI: 1.28 × 10-5-2.73 × 10-5, P value: 6.98 × 10-8), while cg05041783, annotated to gene MARK2, was positively associated with docosapentaenoic acid (DPA) in our fully adjusted model (beta: 9.80 × 10-5, 95%CI: 6.25 × 10-5-1.33 × 10-4, P value: 6.75 × 10-8). In the meta-analysis, we identified the CpG site (cg15951061), annotated to gene CDCA7L below Bonferroni correction (1.23 × 10-7) associated with eicosapentaenoic acid (EPA) intake in model 1 (beta: 2.00 × 10-5, 95% CI: 1.27 × 10-5-2.73 × 10-5, P value = 5.99 × 10-8) and we confirmed the association of cg19937480 with DHA in both models 1 and 2 (beta: 2.07 × 10-5, 95% CI: 1.31 × 10-5-2.83 × 10-5, P value = 1.00 × 10-7 and beta: 2.19 × 10-5, 95% CI: 1.41 × 10-5-2.97 × 10-5, P value = 5.91 × 10-8 respectively). CONCLUSIONS: Our study identified three CpG sites associated with PUFA intake. The mechanisms of these sites remain largely unexplored, highlighting the novelty of our findings. Further research is essential to understand the links between CpG site methylation and PUFA outcomes.

Phenotypic and genome-wide studies on dicarbonyls: major associations to glomerular filtration rate and gamma-glutamyltransferase activity.

BACKGROUND: The dicarbonyl compounds methylglyoxal (MG), glyoxal (GO) and 3-deoxyglucosone (3-DG) have been linked to various diseases. However, disease-independent phenotypic and genotypic association studies with phenome-wide and genome-wide reach, respectively, have not been provided. METHODS: MG, GO and 3-DG were measured by LC-MS in 1304 serum samples of two populations (KORA, n = 482; BiDirect, n = 822) and assessed for associations with genome-wide SNPs (GWAS) and with phenome-wide traits. Redundancy analysis (RDA) was used to identify major independent trait associations. FINDINGS: Mutual correlations of dicarbonyls were highly significant, being stronger between MG and GO (ρ = 0.6) than between 3-DG and MG or GO (ρ = 0.4). Significant phenotypic results included associations of all dicarbonyls with sex, waist-to-hip ratio, glomerular filtration rate (GFR), gamma-glutamyltransferase (GGT), and hypertension, of MG and GO with age and C-reactive protein, of GO and 3-DG with glucose and antidiabetics, of MG with contraceptives, of GO with ferritin, and of 3-DG with smoking. RDA revealed GFR, GGT and, in case of 3-DG, glucose as major contributors to dicarbonyl variance. GWAS did not identify genome-wide significant loci. SNPs previously associated with glyoxalase activity did not reach nominal significance. When multiple testing was restricted to the lead SNPs of GWASs on the traits selected by RDA, 3-DG was found to be associated (p = 2.3 × 10-5) with rs1741177, an eQTL of NF-κB inhibitor NFKBIA. INTERPRETATION: This large-scale, population-based study has identified numerous associations, with GFR and GGT being of pivotal importance, providing unbiased perspectives on dicarbonyls beyond the current state. FUNDING: Deutsche Forschungsgemeinschaft, Helmholtz Munich, German Centre for Cardiovascular Research (DZHK), German Federal Ministry of Research and Education (BMBF).

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Penetrance, cancer incidence and survival in HFE haemochromatosis-A population-based cohort study.

BACKGROUND AND AIMS: Haemochromatosis is characterized by progressive iron overload affecting the liver and can cause cirrhosis and hepatocellular carcinoma. Most haemochromatosis patients are homozygous for p.C282Y in HFE, but only a minority of individuals with this genotype will develop the disease. The aim was to assess the penetrance of iron overload, fibrosis, hepatocellular carcinoma and life expectancy. METHODS: A total of 8839 individuals from the Austrian region of Tyrol were genotyped for the p.C282Y variant between 1997 and 2021. Demographic, laboratory parameters and causes of death were assessed from health records. Penetrance, survival, and cancer incidence were ascertained from diagnosed cases, insurance- and cancer registry data. Outcomes were compared with a propensity score-matched control population. RESULTS: Median age at diagnosis in 542 p.C282Y homozygous individuals was 47.8 years (64% male). At genotyping, the prevalence of iron overload was 55%. The cumulative penetrance of haemochromatosis defined as the presence of provisional iron overload was 24.2% in males and 10.5% in females aged 60 years or younger. Among p.C282Y homozygotes of the same ages, the cumulative proportion of individuals without fibrosis (FIB-4 score < 1.3) was 92.8% in males and 96.7% in females. Median life expectancy was reduced by 6.8 years in individuals homozygous for p.C282Y when compared with population-matched controls (p = .001). Hepatocellular carcinoma incidence was not significantly higher in p.C282Y homozygotes than in controls matched for age and sex. CONCLUSION: Reduced survival and the observed age-dependent increase in penetrance among p.C282Y homozygotes call for earlier diagnosis of haemochromatosis to prevent complications.

Systemic Markers of Lung Function and Forced Expiratory Volume in 1 Second Decline across Diverse Cohorts.

Rationale: Chronic obstructive pulmonary disease (COPD) is a complex disease characterized by airway obstruction and accelerated lung function decline. Our understanding of systemic protein biomarkers associated with COPD remains incomplete. Objectives: To determine what proteins and pathways are associated with impaired pulmonary function in a diverse population. Methods: We studied 6,722 participants across six cohort studies with both aptamer-based proteomic and spirometry data (4,566 predominantly White participants in a discovery analysis and 2,156 African American cohort participants in a validation). In linear regression models, we examined protein associations with baseline forced expiratory volume in 1 second (FEV1) and FEV1/forced vital capacity (FVC). In linear mixed effects models, we investigated the associations of baseline protein levels with rate of FEV1 decline (ml/yr) in 2,777 participants with up to 7 years of follow-up spirometry. Results: We identified 254 proteins associated with FEV1 in our discovery analyses, with 80 proteins validated in the Jackson Heart Study. Novel validated protein associations include kallistatin serine protease inhibitor, growth differentiation factor 2, and tumor necrosis factor-like weak inducer of apoptosis (discovery ß = 0.0561, Q = 4.05 × 10-10; ß = 0.0421, Q = 1.12 × 10-3; and ß = 0.0358, Q = 1.67 × 10-3, respectively). In longitudinal analyses within cohorts with follow-up spirometry, we identified 15 proteins associated with FEV1 decline (Q < 0.05), including elafin leukocyte elastase inhibitor and mucin-associated TFF2 (trefoil factor 2; ß = -4.3 ml/yr, Q = 0.049; ß = -6.1 ml/yr, Q = 0.032, respectively). Pathways and processes highlighted by our study include aberrant extracellular matrix remodeling, enhanced innate immune response, dysregulation of angiogenesis, and coagulation. Conclusions: In this study, we identify and validate novel biomarkers and pathways associated with lung function traits in a racially diverse population. In addition, we identify novel protein markers associated with FEV1 decline. Several protein findings are supported by previously reported genetic signals, highlighting the plausibility of certain biologic pathways. These novel proteins might represent markers for risk stratification, as well as novel molecular targets for treatment of COPD.

Polygenic scores for estimated glomerular filtration rate in a population of general adults and elderly - comparative results from the KORA and AugUR study.

BACKGROUND: Polygenic scores (PGSs) combining genetic variants found to be associated with creatinine-based estimated glomerular filtration rate (eGFRcrea) have been applied in various study populations with different age ranges. This has shown that PGS explain less eGFRcrea variance in the elderly. Our aim was to understand how differences in eGFR variance and the percentage explained by PGS varies between population of general adults and elderly. RESULTS: We derived a PGS for cystatin-based eGFR (eGFRcys) from published genome-wide association studies. We used the 634 variants known for eGFRcrea and the 204 variants identified for eGFRcys to calculate the PGS in two comparable studies capturing a general adult and an elderly population, KORA S4 (n = 2,900; age 24-69 years) and AugUR (n = 2,272, age ≥ 70 years). To identify potential factors determining age-dependent differences on the PGS-explained variance, we evaluated the PGS variance, the eGFR variance, and the beta estimates of PGS association on eGFR. Specifically, we compared frequencies of eGFR-lowering alleles between general adult and elderly individuals and analyzed the influence of comorbidities and medication intake. The PGS for eGFRcrea explained almost twice as much (R2 = 9.6%) of age-/sex adjusted eGFR variance in the general adults compared to the elderly (4.6%). This difference was less pronounced for the PGS for eGFRcys (4.7% or 3.6%, respectively). The beta-estimate of the PGS on eGFRcrea was higher in the general adults compared to the elderly, but similar for the PGS on eGFRcys. The eGFR variance in the elderly was reduced by accounting for comorbidities and medication intake, but this did not explain the difference in R2-values. Allele frequencies between general adult and elderly individuals showed no significant differences except for one variant near APOE (rs429358). We found no enrichment of eGFR-protective alleles in the elderly compared to general adults. CONCLUSIONS: We concluded that the difference in explained variance by PGS was due to the higher age- and sex-adjusted eGFR variance in the elderly and, for eGFRcrea, also by a lower PGS association beta-estimate. Our results provide little evidence for survival or selection bias.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk.

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Identification of risk loci for primary aldosteronism in genome-wide association studies.

Primary aldosteronism affects up to 10% of hypertensive patients and is responsible for treatment resistance and increased cardiovascular risk. Here we perform a genome-wide association study in a discovery cohort of 562 cases and 950 controls and identify three main loci on chromosomes 1, 13 and X; associations on chromosome 1 and 13 are replicated in a second cohort and confirmed by a meta-analysis involving 1162 cases and 3296 controls. The association on chromosome 13 is specific to men and stronger in bilateral adrenal hyperplasia than aldosterone producing adenoma. Candidate genes located within the two loci, CASZ1 and RXFP2, are expressed in human and mouse adrenals in different cell clusters. Their overexpression in adrenocortical cells suppresses mineralocorticoid output under basal and stimulated conditions, without affecting cortisol biosynthesis. Our study identifies the first risk loci for primary aldosteronism and highlights new mechanisms for the development of aldosterone excess.

Uncovering the Contribution of Moderate-Penetrance Susceptibility Genes to Breast Cancer by Whole-Exome Sequencing and Targeted Enrichment Sequencing of Candidate Genes in Women of European Ancestry.

Rare variants in at least 10 genes, including BRCA1, BRCA2, PALB2, ATM, and CHEK2, are associated with increased risk of breast cancer; however, these variants, in combination with common variants identified through genome-wide association studies, explain only a fraction of the familial aggregation of the disease. To identify further susceptibility genes, we performed a two-stage whole-exome sequencing study. In the discovery stage, samples from 1528 breast cancer cases enriched for breast cancer susceptibility and 3733 geographically matched unaffected controls were sequenced. Using five different filtering and gene prioritization strategies, 198 genes were selected for further validation. These genes, and a panel of 32 known or suspected breast cancer susceptibility genes, were assessed in a validation set of 6211 cases and 6019 controls for their association with risk of breast cancer overall, and by estrogen receptor (ER) disease subtypes, using gene burden tests applied to loss-of-function and rare missense variants. Twenty genes showed nominal evidence of association (p-value < 0.05) with either overall or subtype-specific breast cancer. Our study had the statistical power to detect susceptibility genes with effect sizes similar to ATM, CHEK2, and PALB2, however, it was underpowered to identify genes in which susceptibility variants are rarer or confer smaller effect sizes. Larger sample sizes would be required in order to identify such genes.

Integrative analysis of clinical and epigenetic biomarkers of mortality.

DNA methylation (DNAm) has been reported to be associated with many diseases and with mortality. We hypothesized that the integration of DNAm with clinical risk factors would improve mortality prediction. We performed an epigenome-wide association study of whole blood DNAm in relation to mortality in 15 cohorts (n = 15,013). During a mean follow-up of 10 years, there were 4314 deaths from all causes including 1235 cardiovascular disease (CVD) deaths and 868 cancer deaths. Ancestry-stratified meta-analysis of all-cause mortality identified 163 CpGs in European ancestry (EA) and 17 in African ancestry (AA) participants at p < 1 × 10-7 , of which 41 (EA) and 16 (AA) were also associated with CVD death, and 15 (EA) and 9 (AA) with cancer death. We built DNAm-based prediction models for all-cause mortality that predicted mortality risk after adjusting for clinical risk factors. The mortality prediction model trained by integrating DNAm with clinical risk factors showed an improvement in prediction of cancer death with 5% increase in the C-index in a replication cohort, compared with the model including clinical risk factors alone. Mendelian randomization identified 15 putatively causal CpGs in relation to longevity, CVD, or cancer risk. For example, cg06885782 (in KCNQ4) was positively associated with risk for prostate cancer (Beta = 1.2, PMR  = 4.1 × 10-4 ) and negatively associated with longevity (Beta = -1.9, PMR  = 0.02). Pathway analysis revealed that genes associated with mortality-related CpGs are enriched for immune- and cancer-related pathways. We identified replicable DNAm signatures of mortality and demonstrated the potential utility of CpGs as informative biomarkers for prediction of mortality risk.

Colocalization analysis of pancreas eQTLs with risk loci from alcoholic and novel non-alcoholic chronic pancreatitis GWAS suggests potential disease causing mechanisms.

BACKGROUND: Previous genome-wide association studies (GWAS) identified genome-wide significant risk loci in chronic pancreatitis and investigated underlying disease causing mechanisms by simple overlaps with expression quantitative trait loci (eQTLs), a procedure which may often result in false positive conclusions. METHODS: We conducted a GWAS in 584 non-alcoholic chronic pancreatitis (NACP) patients and 6040 healthy controls. Next, we applied Bayesian colocalization analysis of identified genome-wide significant risk loci from both, our recently published alcoholic chronic pancreatitis (ACP) and the novel NACP dataset, with pancreas eQTLs from the GTEx V8 European cohort to prioritize candidate causal genes and extracted credible sets of shared causal variants. RESULTS: Variants at the CTRC (p = 1.22 × 10-21) and SPINK1 (p = 6.59 × 10-47) risk loci reached genome-wide significance in NACP. CTRC risk variants colocalized with CTRC eQTLs in ACP (PP4 = 0.99, PP4/PP3 = 95.51) and NACP (PP4 = 0.99, PP4/PP3 = 95.46). For both diseases, the 95% credible set of shared causal variants consisted of rs497078 and rs545634. CLDN2-MORC4 risk variants colocalized with CLDN2 eQTLs in ACP (PP4 = 0.98, PP4/PP3 = 42.20) and NACP (PP4 = 0.67, PP4/PP3 = 7.18), probably driven by the shared causal variant rs12688220. CONCLUSIONS: A shared causal CTRC risk variant might unfold its pathogenic effect in ACP and NACP by reducing CTRC expression, while the CLDN2-MORC4 shared causal variant rs12688220 may modify ACP and NACP risk by increasing CLDN2 expression.

Epigenetic scores for the circulating proteome as tools for disease prediction.

Protein biomarkers have been identified across many age-related morbidities. However, characterising epigenetic influences could further inform disease predictions. Here, we leverage epigenome-wide data to study links between the DNA methylation (DNAm) signatures of the circulating proteome and incident diseases. Using data from four cohorts, we trained and tested epigenetic scores (EpiScores) for 953 plasma proteins, identifying 109 scores that explained between 1% and 58% of the variance in protein levels after adjusting for known protein quantitative trait loci (pQTL) genetic effects. By projecting these EpiScores into an independent sample (Generation Scotland; n = 9537) and relating them to incident morbidities over a follow-up of 14 years, we uncovered 137 EpiScore-disease associations. These associations were largely independent of immune cell proportions, common lifestyle and health factors, and biological aging. Notably, we found that our diabetes-associated EpiScores highlighted previous top biomarker associations from proteome-wide assessments of diabetes. These EpiScores for protein levels can therefore be a valuable resource for disease prediction and risk stratification.

Although our genetic code does not change throughout our lives, our genes can be turned on and off as a result of epigenetics. Epigenetics can track how the environment and even certain behaviors add or remove small chemical markers to the DNA that makes up the genome. The type and location of these markers may affect whether genes are active or silent, this is, whether the protein coded for by that gene is being produced or not. One common epigenetic marker is known as DNA methylation. DNA methylation has been linked to the levels of a range of proteins in our cells and the risk people have of developing chronic diseases. Blood samples can be used to determine the epigenetic markers a person has on their genome and to study the abundance of many proteins. Gadd, Hillary, McCartney, Zaghlool et al. studied the relationships between DNA methylation and the abundance of 953 different proteins in blood samples from individuals in the German KORA cohort and the Scottish Lothian Birth Cohort 1936. They then used machine learning to analyze the relationship between epigenetic markers found in people's blood and the abundance of proteins, obtaining epigenetic scores or 'EpiScores' for each protein. They found 109 proteins for which DNA methylation patterns explained between at least 1% and up to 58% of the variation in protein levels. Integrating the 'EpiScores' with 14 years of medical records for more than 9000 individuals from the Generation Scotland study revealed 130 connections between EpiScores for proteins and a future diagnosis of common adverse health outcomes. These included diabetes, stroke, depression, various cancers, and inflammatory conditions such as rheumatoid arthritis and inflammatory bowel disease. Age-related chronic diseases are a growing issue worldwide and place pressure on healthcare systems. They also severely reduce quality of life for individuals over many years. This work shows how epigenetic scores based on protein levels in the blood could predict a person's risk of several of these diseases. In the case of type 2 diabetes, the EpiScore results replicated previous research linking protein levels in the blood to future diagnosis of diabetes. Protein EpiScores could therefore allow researchers to identify people with the highest risk of disease, making it possible to intervene early and prevent these people from developing chronic conditions as they age.

Genome-wide association study identifies susceptibility loci for acute myeloid leukemia.

Acute myeloid leukemia (AML) is a hematological malignancy with an undefined heritable risk. Here we perform a meta-analysis of three genome-wide association studies, with replication in a fourth study, incorporating a total of 4018 AML cases and 10488 controls. We identify a genome-wide significant risk locus for AML at 11q13.2 (rs4930561; P = 2.15 × 10-8; KMT5B). We also identify a genome-wide significant risk locus for the cytogenetically normal AML sub-group (N = 1287) at 6p21.32 (rs3916765; P = 1.51 × 10-10; HLA). Our results inform on AML etiology and identify putative functional genes operating in histone methylation (KMT5B) and immune function (HLA).

Rare and low-frequency exonic variants and gene-by-smoking interactions in pulmonary function.

Genome-wide association studies have identified numerous common genetic variants associated with spirometric measures of pulmonary function, including forced expiratory volume in one second (FEV1), forced vital capacity, and their ratio. However, variants with lower minor allele frequencies are less explored. We conducted a large-scale gene-smoking interaction meta-analysis on exonic rare and low-frequency variants involving 44,429 individuals of European ancestry in the discovery stage and sought replication in the UK BiLEVE study with 45,133 European ancestry samples and UK Biobank study with 59,478 samples. We leveraged data on cigarette smoking, the major environmental risk factor for reduced lung function, by testing gene-by-smoking interaction effects only and simultaneously testing the genetic main effects and interaction effects. The most statistically significant signal that replicated was a previously reported low-frequency signal in GPR126, distinct from common variant associations in this gene. Although only nominal replication was obtained for a top rare variant signal rs142935352 in one of the two studies, interaction and joint tests for current smoking and PDE3B were significantly associated with FEV1. This study investigates the utility of assessing gene-by-smoking interactions and underscores their effects on potential pulmonary function.

Plasma Proteomics of Renal Function: A Transethnic Meta-Analysis and Mendelian Randomization Study.

BACKGROUND: Studies on the relationship between renal function and the human plasma proteome have identified several potential biomarkers. However, investigations have been conducted largely in European populations, and causality of the associations between plasma proteins and kidney function has never been addressed. METHODS: A cross-sectional study of 993 plasma proteins among 2882 participants in four studies of European and admixed ancestries (KORA, INTERVAL, HUNT, QMDiab) identified transethnic associations between eGFR/CKD and proteomic biomarkers. For the replicated associations, two-sample bidirectional Mendelian randomization (MR) was used to investigate potential causal relationships. Publicly available datasets and transcriptomic data from independent studies were used to examine the association between gene expression in kidney tissue and eGFR. RESULTS: In total, 57 plasma proteins were associated with eGFR, including one novel protein. Of these, 23 were additionally associated with CKD. The strongest inferred causal effect was the positive effect of eGFR on testican-2, in line with the known biological role of this protein and the expression of its protein-coding gene (SPOCK2) in renal tissue. We also observed suggestive evidence of an effect of melanoma inhibitory activity (MIA), carbonic anhydrase III, and cystatin-M on eGFR. CONCLUSIONS: In a discovery-replication setting, we identified 57 proteins transethnically associated with eGFR. The revealed causal relationships are an important stepping stone in establishing testican-2 as a clinically relevant physiological marker of kidney disease progression, and point to additional proteins warranting further investigation.

Metabolic syndrome and the plasma proteome: from association to causation.

BACKGROUND: The metabolic syndrome (MetS), defined by the simultaneous clustering of cardio-metabolic risk factors, is a significant worldwide public health burden with an estimated 25% prevalence worldwide. The pathogenesis of MetS is not entirely clear and the use of molecular level data could help uncover common pathogenic pathways behind the observed clustering. METHODS: Using a highly multiplexed aptamer-based affinity proteomics platform, we examined associations between plasma proteins and prevalent and incident MetS in the KORA cohort (n = 998) and replicated our results for prevalent MetS in the HUNT3 study (n = 923). We applied logistic regression models adjusted for age, sex, smoking status, and physical activity. We used the bootstrap ranking algorithm of least absolute shrinkage and selection operator (LASSO) to select a predictive model from the incident MetS associated proteins and used area under the curve (AUC) to assess its performance. Finally, we investigated the causal effect of the replicated proteins on MetS using two-sample Mendelian randomization. RESULTS: Prevalent MetS was associated with 116 proteins, of which 53 replicated in HUNT. These included previously reported proteins like leptin, and new proteins like NTR domain-containing protein 2 and endoplasmic reticulum protein 29. Incident MetS was associated with 14 proteins in KORA, of which 13 overlap the prevalent MetS associated proteins with soluble advanced glycosylation end product-specific receptor (sRAGE) being unique to incident MetS. The LASSO selected an eight-protein predictive model with an (AUC = 0.75; 95% CI = 0.71-0.79) in KORA. Mendelian randomization suggested causal effects of three proteins on MetS, namely apolipoprotein E2 (APOE2) (Wald-Ratio = - 0.12, Wald-p = 3.63e-13), apolipoprotein B (APOB) (Wald-Ratio = - 0.09, Wald-p = 2.54e-04) and proto-oncogene tyrosine-protein kinase receptor (RET) (Wald-Ratio = 0.10, Wald-p = 5.40e-04). CONCLUSIONS: Our findings offer new insights into the plasma proteome underlying MetS and identify new protein associations. We reveal possible casual effects of APOE2, APOB and RET on MetS. Our results highlight protein candidates that could potentially serve as targets for prevention and therapy.

Epigenome-wide association study of whole blood gene expression in Framingham Heart Study participants provides molecular insight into the potential role of CHRNA5 in cigarette smoking-related lung diseases.

BACKGROUND: DNA methylation is a key epigenetic modification that can directly affect gene regulation. DNA methylation is highly influenced by environmental factors such as cigarette smoking, which is causally related to chronic obstructive pulmonary disease (COPD) and lung cancer. To date, there have been few large-scale, combined analyses of DNA methylation and gene expression and their interrelations with lung diseases. RESULTS: We performed an epigenome-wide association study of whole blood gene expression in ~ 6000 individuals from four cohorts. We discovered and replicated numerous CpGs associated with the expression of cis genes within 500 kb of each CpG, with 148 to 1,741 cis CpG-transcript pairs identified across cohorts. We found that the closer a CpG resided to a transcription start site, the larger its effect size, and that 36% of cis CpG-transcript pairs share the same causal genetic variant. Mendelian randomization analyses revealed that hypomethylation and lower expression of CHRNA5, which encodes a smoking-related nicotinic receptor, are causally linked to increased risk of COPD and lung cancer. This putatively causal relationship was further validated in lung tissue data. CONCLUSIONS: Our results provide a large and comprehensive association study of whole blood DNA methylation with gene expression. Expression platform differences rather than population differences are critical to the replication of cis CpG-transcript pairs. The low reproducibility of trans CpG-transcript pairs suggests that DNA methylation regulates nearby rather than remote gene expression. The putatively causal roles of methylation and expression of CHRNA5 in relation to COPD and lung cancer provide evidence for a mechanistic link between patterns of smoking-related epigenetic variation and lung diseases, and highlight potential therapeutic targets for lung diseases and smoking cessation.

Genome-wide association study identifies novel susceptibility loci for KIT D816V positive mastocytosis.

Mastocytosis is a rare myeloid neoplasm characterized by uncontrolled expansion of mast cells, driven in >80% of affected individuals by acquisition of the KIT D816V mutation. To explore the hypothesis that inherited variation predisposes to mastocytosis, we performed a two-stage genome-wide association study, analyzing 1,035 individuals with KIT D816V positive disease and 17,960 healthy control individuals from five European populations. After quality control, we tested 592,007 SNPs at stage 1 and 75 SNPs at stage 2 for association by using logistic regression and performed a fixed effects meta-analysis to combine evidence across the two stages. From the meta-analysis, we identified three intergenic SNPs associated with mastocytosis that achieved genome-wide significance without heterogeneity between cohorts: rs4616402 (pmeta = 1.37 × 10-15, OR = 1.52), rs4662380 (pmeta = 2.11 × 10-12, OR = 1.46), and rs13077541 (pmeta = 2.10 × 10-9, OR = 1.33). Expression quantitative trait analyses demonstrated that rs4616402 is associated with the expression of CEBPA (peQTL = 2.3 × 10-14), a gene encoding a transcription factor known to play a critical role in myelopoiesis. The role of the other two SNPs is less clear: rs4662380 is associated with expression of the long non-coding RNA gene TEX41 (peQTL = 2.55 × 10-11), whereas rs13077541 is associated with the expression of TBL1XR1, which encodes transducin (ß)-like 1 X-linked receptor 1 (peQTL = 5.70 × 10-8). In individuals with available data and non-advanced disease, rs4616402 was associated with age at presentation (p = 0.009; beta = 4.41; n = 422). Additional focused analysis identified suggestive associations between mastocytosis and genetic variation at TERT, TPSAB1/TPSB2, and IL13. These findings demonstrate that multiple germline variants predispose to KIT D816V positive mastocytosis and provide novel avenues for functional investigation.

Meta-analysis uncovers genome-wide significant variants for rapid kidney function decline.

Rapid decline of glomerular filtration rate estimated from creatinine (eGFRcrea) is associated with severe clinical endpoints. In contrast to cross-sectionally assessed eGFRcrea, the genetic basis for rapid eGFRcrea decline is largely unknown. To help define this, we meta-analyzed 42 genome-wide association studies from the Chronic Kidney Diseases Genetics Consortium and United Kingdom Biobank to identify genetic loci for rapid eGFRcrea decline. Two definitions of eGFRcrea decline were used: 3 mL/min/1.73m2/year or more ("Rapid3"; encompassing 34,874 cases, 107,090 controls) and eGFRcrea decline 25% or more and eGFRcrea under 60 mL/min/1.73m2 at follow-up among those with eGFRcrea 60 mL/min/1.73m2 or more at baseline ("CKDi25"; encompassing 19,901 cases, 175,244 controls). Seven independent variants were identified across six loci for Rapid3 and/or CKDi25: consisting of five variants at four loci with genome-wide significance (near UMOD-PDILT (2), PRKAG2, WDR72, OR2S2) and two variants among 265 known eGFRcrea variants (near GATM, LARP4B). All these loci were novel for Rapid3 and/or CKDi25 and our bioinformatic follow-up prioritized variants and genes underneath these loci. The OR2S2 locus is novel for any eGFRcrea trait including interesting candidates. For the five genome-wide significant lead variants, we found supporting effects for annual change in blood urea nitrogen or cystatin-based eGFR, but not for GATM or LARP4B. Individuals at high compared to those at low genetic risk (8-14 vs. 0-5 adverse alleles) had a 1.20-fold increased risk of acute kidney injury (95% confidence interval 1.08-1.33). Thus, our identified loci for rapid kidney function decline may help prioritize therapeutic targets and identify mechanisms and individuals at risk for sustained deterioration of kidney function.